Observation of a mixed pairing state in Ru microdomains embedded in Sr2Ru04 
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We report first measurements on tunneling into the interior of a Ru microdomain embedded in 
bulk chiral p-wave superconductor S^RuCU. The junctions were prepared by pressing a pure In 
wire onto the ab face of a Ru-containing S^RuCU single crystal. Below the superconducting tran- 
sition temperature (T c ) of S^RuCU (= 1.5 K), we observed a zero-bias conductance peak (ZBCP) 
associated with the proximity-induced p-wave superconductivity. However, below 0.44 K, roughly 
the T c of bulk Ru (= 0.5 K), an unexpected Ru gap and a broad ZBCP superimposed with subgap 
features were found to emerge, with unusual magnetic field dependences. We argue that the new 
features resulted from a mixed pairing state featuring intrinsic s- and p-wave condensates in the 
interior of a Ru microdomain embedded in Sr2Ru04. 

PACS numbers: 74.50.+r, 74.25.Fy, 74.70.-b 



Superconductors are divided into two categories fea- 
turing an even-parity, spin-singlet and odd-parity, spin- 
triplet pairing symmetry, respectively. The inversion 
symmetry requires that the pairing symmetry of a bulk 
superconductor be of either an even or odd parity. How- 
ever, in the absence of such symmetry, as in crystals with 
no inversion center or near a superconductor-normal 
metal interface 0, mixed pairing states are allowed by 
symmetry, and realizable under suitable conditions. 

Pairing symmetry of a superconductor is determined 
primarily by effective interactions. In the presence of 
a rotational symmetry, the interaction can be decom- 
posed into channels with different orbital angular mo- 
mentum, I. The superconducting energy gap, A;, is given 
by A/ = 2eiexp(— 2/NqVi), where £j is a characteristic 
energy, No is the density of states (DOS) at the Fermi 
energy, and Vi is interaction in the I channel 3] . States of 
I = 0, 2 ... then correspond to s-, rf-wave ...(even-parity) 
and those of I = 1, 3 ... to p-, /-wave ... (odd-parity) pair- 
ings. Electrons will usually pair in the channel with the 
largest Vi , corresponding to the highest superconducting 
transition temperature (T c ), to attain a largest conden- 
sation energy, even though other V/'s may be substantial. 

Chiral p-wave superconductivity has been established 
in Si'2Ru04 by many experiments 0] including recent 
phase sensitive measurements The eutectic Ru- 

Sr2Ru04 system of pure Ru microdomains (a few ^m 
in size) embedded in bulk Sr2Ru04 |fj|], which is of fun- 
damental interest because of the enhanced T c as high 
as 3 K, has emerged as a useful system for the study of 
unconventional pairings 0, . Bulk Ru is an s-wave su- 
perconductor, as evidenced by the insensitivity of super- 
conductivity to impurities |l0j. The intriguing questions 
are how the intrinsic pairing in s-wave channel competes 
for stability with the proximity-induced p-wave pairing 
below the intrinsic T c of Ru and whether intrinsic pair- 



ing in p-wave channel is possible in the interior of a Ru 
microdomain embedded in Sr2Ru04. 

To answer these questions, we prepared junctions of 
In/Ru/Sr2Ru04 on S^RuCm single crystals grown by 
floating-zone technique 4]. These crystals are essen- 
tially pure crystals containing only few scattered Ru mi- 
crodomains. A Sr2RuC>4 crystal was first cleaved along 
its ab plane. Then a freshly cut, high-purity (5N) In 
wire of a 0.25 mm diameter was pressed onto the ab 
face of the crystal without any artificial tunnel barrier. 
The current-voltage (I-V) characteristics of the junctions 
were measured with a standard 4-terminal, DC technique 
in a 3 He and dilution refrigerator with a base tempera- 
ture of 300 and 20 mK, respectively. The tunneling con- 
ductance di~/dV were obtained by taking the numeric 
derivatives of the I-V curve. 

Even though no structural characterization of the junc- 
tions was carried out, the totality of the tunneling and 
transport data allows us to divide the junctions into three 
types. For the first type, the tunneling is dominated by a 
contact between In and the nonsuperconducting surface 
layer 0] of pure Sr 2 Ru04, yielding a tunnel spectrum 
well described by the Bardeen-Cooper-Schiffer (BCS) for- 
mula 01 with an In gap of 0.5 meV. For the second 
type, the tunneling is dominated by a contact between 
In and the 3-K phase, the interface region between Ru 
and Sr2Ru04, showing an In gap around 0.5 meV and 
a zero-bias conductance peak (ZBCP) resulted from An- 
dreev surface bound states (ASBSs), similar to those seen 
previously |8|,|9(- For the third type of junctions, the fo- 
cus of the present work, the tunneling is dominated by a 
contact between In and the interior of a Ru microdomain 
(Inset of Fig. la). 

Measurements on the temperature-dependence of 
the junction resistance, Rj(T), revealed a negative 
dRj(T)/dT above 3.4 K (Fig. la), suggesting the ex- 
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FIG. 1: (color online) (a) Rj(T)/Rj(5K) vs T for two 
In/Ru/Sr2Ru04 junctions. InRu-1 was measured in three 
separate cooling-downs - in a 3 He refrigerator with a c- 
axis field (InRu-lA) and a dilution refrigerator with a c- 
axis (InRu-lB) and an in-plane (InRu-lC) field. Values for 
Rj{5K) are 47.5, 51.2 and 60.6 fi for InRu-lA, B, and C 
and 12.5 SI for InRu-2. Inset: Schematic of the third type of 
In/Ru/Sr2RuC>4 junctions. The nonsuperconducting surface 
layer (N) is indicated; (b) dl/dV vs V for the InRu-lC for 
T > 0.5 K in the zero field. Curves except the bottom one 
(1.6 K) are shifted upward by 0.01 for clarity. 



istence of a barrier between In and Ru, probably due 
to the formation of InOa; on In surface, resulting in 
a superconductor-insulator-superconductor (SIS') junc- 
tion. Rj(T) showss a drop at the T c of In, 3.4 K, for 
all samples, and a barely visible kink around 1.5 K, the 
T c of bulk Si'2Ru04 for all samples except InRu-lC. For 
the latter, a much more pronounced kink was found at 
a slightly higher temperature (1.7 K). We believe that 
the shift is due to the development of additional contacts 
between In and Ru close to, but not at the 3-K phase 
region after two thermal cyclings. No feature was found 



between 1.7 to 3 K, suggesting that the In contact is away 
from the 3-K phase in all samples. As the temperature 
is lowered, Rj (T) exhibits an abrupt drop around the T c 
of bulk Ru, 0.5 K. For InRu-lB and C, a non-monotonic 
Rj(T), to be discussed below, is seen below 0.5 K. 

Tunnel spectra T > 0.5 K are shown in Fig. lb for 
InRu-lC. Similar spectra were obtained for other sam- 
ples. A rather small ZBCP was found as high as 1.6 K, 
above the T c of bulk Sr2Ru04. Such ZBCP is very similar 
to those found previously in tunneling into the 3-K phase 
region where Andreev surface bound states (ASBS) are 
formed on the sample surface of the 3-K phase • These 
ASBSs are formed in this case because of an exponen- 
tially small, proximity-induced p-wave order parameter 
in the interior of the Ru microdomain. As the sample 
was cooled to below 0.44 K, slightly lower than the T c 
of bulk Ru (0.5 K), new features emerge in the tunnel 
spectrum, including a sharp drop in the tunneling con- 
ductance at a well-define gap edge, Ar u (=qVr u ), and 
a ZBCP (Fig. 2a). Again almost identical features were 
found for all other samples. 

For such SIS' junction, gap features would normally 
be expected at A/ n ± Ar u . However, the spectrum can 
change when broadening is introduced. Starting from the 
normalized I — V characteristic for a SIS' junction, 



I(V) 



N(E)N'(E+V)[f(E)-f(E+eV)]dE (1) 



where E is the energy, f(E) is the Fermi function, 
N(E) and N'(E) are the normalized, broadened DOS 
for S and S' of the form 



N(E) = Re[ 



E I +iT 



(2) 



where T is the broadening parameter reflecting, for ex- 
ample, a finite quasiparticle lifetime |l3|. we can obtain 
the normalized tunnel conductance Gj (= dl/dV). As- 
suming that In and Ru are s-wave superconductors, we 
indeed obtain gap features at A#„/e and Aj n /e in our 
simulation, as shown in Fig. 2b. Interes tingly, similar 
results were obtained for high-T c junctions 141.1 15|. Iden- 
tifying Ar u as the Ru gap is further supported by the 
excellent fit of A Ru (T) to the BCS formula (Fig. 2c). 

Another important observation is that the ZBCP 
evolves into a continuous plateau superimposed by dis- 
tinct peaks below 0.44 K. The conductance values of 
the plateau were found to be larger than twice of the 
normal-state conductance, characteristic of unconven- 
tional pairing, suggesting that p-wave pairing continues 
to be present at low temperatures. The satellite peak, 
found around Ar u /2 in the ZBCP just below 0.44 K 
and seen to split into two at the lowest temperatures, 
should correspond to subgap resonant states. The sharp 
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FIG. 2: (a) dJ/dV vs V for the InRu-lC for T < 0.5 K in 
zero field. Curves except the bottom one are shifted upward 
by 0.03 n _1 for clarity. The red arrows indicate Vr u \ (b) 
Normalized AI /dV vs V for the InRu-lC in zero field at 0.1 
K (upper curve with A In ^ Ru ) = 0.50 (0.10) meV). The lower 
curve shows a numerical simulation of a conventional In/Ru 
junction at T = 0.1 K (shifted by 0.6 fi -1 ), employing an 
effective DOS of the form, N eff (E) = r(N(E) -1) + 1), where 
N(E) is the broadened BCS DOS (Eq. 2) and r accounts for 
the junction leakage (N e ff(E) = 1 at high energies). The 
parameters used in the simulation are: A In ^ Ru -) = 0.50 (0.072) 
meV, r In{Ru) = 0.20 (0.03) meV, and r In(Ru) = 0.15 (0.50); 
(c) Un-normalized A_r„ vs T for all samples. Solid lines are 
BCS fits to A Ru vsT, yielding A Ru = 0.097, 0.097, 0.104, and 
0.108 meV, and T c = 0.44, 0.44, 0.44 and 0.45 K for InRu-lA, 
B, C, and InRu-2, respectively. 



central peak in the spectrum, however, must come from 
the Josephson coupling between In and Ru. If the order 
parameter direction of proximity-induced p-wave state is 
the same as that in the bulk, as expected, Josephson 
coupling between the s-wave In and proximity-induced 
p-wave Ru is forbidden by the selection rule [lfj . There- 
fore the observed Josephson coupling may indicate the 
emergence of a new superconducting component. 

The development of a gap edge close to the T c of 
bulk Ru suggests that the intrinsic s-wave pairing of Ru 
may be this new superconducting component, suggesting 
father pairings in both I = and 1 channels in the inte- 
rior of a Ru microdomain below 0.44 K. The interesting 
question is whether the p-wave pairing can be intrinsic, 
rather than proximity-induced. Such simultaneous pair- 
ings in both s- and p-wave channels in a Ru microdomain 
is not forbidden by symmetry. The typical size of a Ru 
microdomain is a few /im. The superconducting coher- 
ence length of bulk Ru, £,o.Ru, in the clean limit, is w 4/zm 
based on (,o,Ru — l^^-o,Ru, where vf is the Fermi ve- 
locity and A 0j fl M is energy gap (vf = 1.5 x 10 6 m/s and 
A.o,r u = 0.072 meV). An inversion operation on such 
a Ru microdomain does not constrain the wavefunction 
symmetry of the Cooper pair as its size is comparable 
with or larger than the sample size. Microscopically, Vi 
usually decreases rapidly with increasing I, making Vi 
the second largest among VJs. A condensate in the I = 1 
channel in the Ru microdomain, with a smaller conden- 
sation energy than that of a I = condensate, could be 
stabilized by energy saving from minimizing the gradient 
of the order parameter between the Ru microdomain and 
the 3-K phase of Sr 2 Ru04. 

This scenario is consistent with the quantitative result 
on Ao.7?u (= 0.1 meV), which is larger than the BCS 
value, 0.072 meV, estimated from T c = 0.44 K. The ra- 
tio, 2A /kBT c , has an unexpectedly high value (w 5.5), 
typical of a strong-coupling superconductor, which is dif- 
ficult to accept because of the low T c . In fact, from 
T c = l.UT g exp[-l/N{0)V] H3 and Tg :Ru w 404 K, 
we find that [^(0)^]^ — 0.147, comparable with that 
of a protypical weak-coupling superconductor, Al, with 
[iV(0)V]Az ~ 0.165. On the other hand, if we assume a 
mixed pairing state with the form |l7j| : 



A 



Ru 



iA r 



(3) 



where A s — ia^Ao is the s-wave and A p = ioi<j ■ d the 
p-wave gap function, with d the vector order parameter 
and <7 = (<7i, era, 03) Pauli matrices |3j. If A p belongs to 
a one-dimensional representation and is real with a mag- 



nitude similar to that of A^ 
we find that 



A 



+ 



\ p |w 0.072 meV, 
F = 0.102 meV, 



which agrees well with experiment. This scenario also 
provides a plausible explanation for the non-monotonic 
behavior in Rj(T), seen also in InRu-2 under a magnetic 
field. Rj(T) of an SIS' junction measured above its crit- 
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FIG. 3: dl/dV vs V for InRu-1 (T = 0.1 K) with the field ap- 
plied along the ab plane (a) and the c axis (b). Curves except 
the bottom ones are shifted upward by 0.035 SI -1 for clarity. 
The red arrows track the field variation of Ar u . For field ap- 
plied along the c axis, the spectrum depends on whether the 
sample is zero-field- or field-cooled. Data presented here were 
obtained after the sample was field-cooled from 5 K. 



ical current, I c , follows the temperature dependence of 
J C (T). The Josephson couplings between In and the two 
individual superconducting components in Ru could be 
of different signs, leading to a non-monotonic I C (T) and 
Rj (T) , as seen previously 0] . 

The magnetic field dependence of tunnel spectrum pro- 
vides further support to the mixed pairing picture. As- 
suming that the field dependence of A s is essentially that 
in the bulk, a type I superconductor with a H c of 69 Oe, 
A s will be fully suppressed at 69 Oe. As shown in Fig. 3, 
with a field applied along the ab plane, Ar u is seen to be 
suppressed rapidly as the field increase to about 70 Oe 
and more slowly for larger field. Within the mixed pair- 
ing picture, A s vanishes at this field, leading to Ar u = 
A p . Experimentally, A^ u = 0.07 meV at 70 Oe, agreeing 
with the mixed pairing picture. 



An intriguing evolution of tunneling spectrum was ob- 
served when the field was applied along the c axis. As 
seen in Fig. 3b, Ar u shows a non-monotonic behavior - 
with the increasing field, the gap and ZBCP features seen 
at lower fields were found to re-emerge at higher fields. 
Qualitatively similar behavior was observed in InRu-lA 
and in InRu-2. Since Sr2Ru04 is a chiral p-wave super- 
conductor [4j, chiral currents flowing in the ab plane are 
expected on the surface of the Ru microdomain. When 
the interior of the Ru microdomain becomes supercon- 
ducting, a chiral order parameter given by Eq. 3 will 
lead to its own chiral current. When a magnetic field is 
applied along the c axis, a finite flux will thread through 
the in-plane chiral current loop. We believe that the flux 
will modulate the chiral current just as the flux modu- 
lates the supercurrent in Little-Parks effect 01 > leading 
to T c and gap modulations. This process will minimize 
(maximize) the gap edge (Ar u ) for certain fields, as seen 
experimentally. Such modulation will not occur when the 
field is applied along the in-plane direction. 
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